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CHAPTER 1 

The Hoyle-Narlikar Theor y 

of Gravitation 


1. Introduction 


The success of Maxwell's equations has led to 
electrodynamics being normally formulated in terms of fields 
that have degrees of freedom independent of the particles in 
them. However, Gauss suggested that an action-at-a-distance 
theory in which the action travelled at a finite velocity 
might be possible. This idea was developed by Wheeler and 
Feynman ^who derived their theory from an action-principl 
that involved only direct interactions between pairs oj part— 
icles. a feature of this theory was that the 'pseudo'-fields 


introduced are the half-retarded plus half-advanced field 






laculated from the world-lines of the particles. However, 

(z) 

Wheeler and Feynman, and, in a different way, Hogarth 
were able to show that, provided certain cosmological 


conditions were satisfied, these fields could combine to 


give the observed field. Hoyle and Karlikar 


( 4 ) 


extended the 


ttieory to general space—times and oooained similar theories 
for their 'C’-field ^^and f 0 r the gravitational fabeld v °'. 


It is with these theories that this chapter is concerned. 






















It will be . 5 'aown that in an expanding universe the 


advanced fields are infinite,and the retarded fields finite. 
This is because, unlike electric charges, all masses have the 
same sign. 


2. The Boundary Condition 


Hoyle and Narlikar derive their theory from the 


action: 


4 



Q d'O. d b. 


A f b 


where the 

b. , ■ . 


a 


integration is over the world-lines of particles 
In this expression ^ is a Green function . 


that satisfies the wave equation: 


C (x, x') j 


* * 


^Q(xX) * fVx. x 


where 3 is the determinant of . Since the double sum 

in the action A is symmetrical between all pairs of 
particles a f b , only that part of Qf<x,b) that is 

symmetrical between cu and ti will contribute to tne action 
i,e. the action can be written 

A " 

where ^ ( 




CK 


^ jj C C<X,b) cta y cl Jb 

b) r 2 «- (OLj b ) V a Q (b* °0. 


m hus ({ must be the time-symmetric Green function, and can 


X 


be written 


<?*. i <? 


ret 




where 


G 


A't? tr 
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and 


are the retarded and advanced Green functions. 


^c-c) J 


6j reouinn 


that the action be stationary under variations 


of the 


[ 




^ b 


Q:i , Hoyle and Karlikar obtain the field-equations: 






(b) \ a / (a) (b) 


fix. 


• i « 


), 2 C 


n ., ox 


; * 


- ‘*5^ 


/vi 


M ; r (b) 


M 


r 


l 


w 


here Ax^(x) - j~Q*Yx, a)la. 


However, as a 


consequence of the particular choice of Green function, the 
contraction of the field-equations is satisfied identically, 
'there are thus only 9 equations for the 10 components of c ( -/ 
and the system is indeterminate. 

iioyle and Narlikar therefore impose X m - =const.. 


as the tenth equation. By then making the 'smooth' 

/ < ( a ) (jo) 

approximation, that is by putting „7-h_/Vv fix 

a £ b 

they obtain the Einstein field-equations: 


fluid' 

2 


fi\ 


O ) 


6 Av 


3 


b- 


) 


- / 


0 K 


fhereis an important difference, however, between these 

field-equations in the direct-particle interaction tceory 

and in the usual general theory of relativity. In the 
general theory of relativity, any metric that satisfies the 
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the field-equations is admissible, but in the direct-particle 
interaction theory only those solutions of the field-equations 
are admissible that satisfy the additional requirement: 






This requirement is highly restrictive; it will be shown 
that it is not satisfied for the cosmological solutions of 
the minstein field-equations, and it appears that it cannot 
be satisfied for any models of the universe that either 

i 

contain an infinite amount of matter or undergo infinite 
expansion. 

The difficulty is similar to that occurring in 
Newtonian theory when it is recognized that the universe 
might be infinite. 

The Newtonian potential (p obeys the equation: 

where p is the density. 




W 
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In an, infinite static universe 
the source always has the same 


(b woulc 


lu be infinite, since 


sign. The difficulty was resol¬ 


ved when it was realized that the universe was expanding, since 
in an expanding universe the retarded solution of the above 
equation is finite by a sort of' red-shift ' effect. The 
advanced solution will be infinite by a 'blue-shift' effect. 
This is unimportant in Newtoni cin thoox’y, since one is free 
to choose the solution of the equation and so may ignore the 
infinite advanced solution and take simply the finite 
retarded solution. 

similarly in the direct-particle interaction theory the 


m - field satisfies the equation: 

On * Z ,JVL r A J 


( A/ ?o) / 


where A/ is the density of world-lines of particles. 


As in 


the Newtonian case, one may expect that the effect of the 
expansion of the universe will be to make the retarded solution 
finite and the advanced solution infinite. However, one is 
now not free to choose the finite retarded solution, for the 
equation is derived from a direct-particle interaction action- 
principle symmetric between pairs of particles, and one must 
choose for M, half the sum of the retarded and advanced 
solutions. We would expect tais to be infinite, and this is 
shown to be so in the next section. 
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b. The Cosmological Solutions 


The Robertson-Walker cosmological metrics have the 


t) 


form 


is h id 


R 2 U) 


I _ 2 

ctr 


/ - f\ y- 2* 


v r i (d.Q Z -* 


La 


Since they are conformally flat, one can choose coordinate 


in which they become 

ds 2 - Jl* Let? 
- J7* 


2 


d 


f 


a 


7 a t 


dx 0, dz 


t 


T P ^ A 


2 


$ ci 



w here 


7 


^ is the flat-space metric tensor and 


hi 


Ifpjfl) - 





( 7 ) 


/ffj +T« ) a J j 


For example, for the Finstein-de Sitter -universe 

K.O d (t) - fb\ 5 ( 0 <. 1 ^ oO ) 



JT * k = /X w 

(' ”2 7 J t 


(© t < oO) 


/ 


^ ’*) 


r - r 

For the steady-state (de Sitter) universe 

K 0 kit) s' ■_t (-co-ct^oo) 

e-T 
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JS 


H 


- i 


( - 


<- c 


<- O 


r 


r 




Te ) 


Q (a, l>) obeys the equation 

-\-*k 


SI 


The Green function 

a q% k)-* t -- CCa,b) 

* + ——_— 

V — Q 

r'roia this it follows that J 

L <) _/ry 1, 1 Q 

r / jv‘ 


4 


o> x 


a 


X 







= Jl 


H $''(a, b) 


If we let Q*.n-S. 


t:ien 


SL ^ 


<3 


X 


(X 



T“B 

c?x 


s 


f V i) 


This is simply the flat-space Green function equation, ano. 


aence 


Si'cpj t, p) -- 




m 


he Vv.-field is given by 

*.(*,) ScQ*A/S- 



clx'* - i ( 


"W r * 


\ 

)„ > 


O-rJi/ 


-?or universes wittiout cree-iioxi (e * g. "bne -Litinstein^de £)iot-8 

AJ lZ ~ 3 


universe), 


rt . 


rv - 


const. -c or 
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universes 


const:. 


with creation (steady state) (V - fi /I 

\ W * / 




) __ Jl" (? •) r A fJKt a) G <rr cU 


C rrr 


where the integration is over the future light cone, ahis 
will normally be infinite in an expanding universe, e.g. in 
the iiinstein-de fitter universe. 






oO 


\ 

_i 
T / 


rt 


c 


t 


't 


) n 


? 






- oO 


In the steady-state universe 


A 

do 1 




<0 


- I 


t 


] 


-i rO 


rt 


J 



( 




) 


&o 


By contrast, on the other hand, we have 


(K. 


ttb 


) 


s J7 'Va) 



NJl Urr'r Z cL^ 

(4 rrt 


where the integration is over the past light cone. This wil 
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Inches 
I 


I I i 


cm 


1 



1 I 2 


normally be finite, e.g. in the i)instein~de Sitter universe 

'2 /*■ t. 


(K 


Kit. 


it,) 


— rv 


C 


- jfU 


while in the steady-state universe 

- f 


‘“O - 


n 


— cO 


-r \ 3 




Z 


z 


z n I 


Thus it can be seen that the solution ~ const, of the 


equation 


— n 

□ <n. uh 


is not, in a cosmological metric, the half-advanced plus 
half-retarded solution since this would be infinite. In fact, 
in. the case of the iiinstein-de Sitter and steady-state metrics, 
it is the pure retarded solution. 


4. The 'G 1 -Fi 


Hoyle and Marlifcar derive their direct-particle 
interaction theory of the 'O'-field from the action 

=■ £ £ ff Q(a, i>)-i k <fc£eU> * 

(X fb ' 


© Stephen Hawking, All Rights Reserved 

















where the suffixes a b refer to differentiation of 


Cj (a , b) 


on the world-lines of CK, b respectively. 


t f is a Green function obeying the equation 

dq (x,X') . P(x,X') 




5 


We define the 'G'-field by 


C(o c) = £ fd . W 



a 


OL 


and the matter-current 


J 


v K 


by 


k 




(j) = S k ( y, b) ci b , 


Then 


C 60 


ccr 


X 


1 


h K ^ 


3 




X 


J 


ac 


j 


r K 


K 


We thus see that the sources of the 'O'—field are the places 

where matter is created or destroyed. 

As in the case of the '/d-' -field, the Green function 

must be time-symmetric, that is 


(J (*,0 




I 
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Hoyle and Narlikar claim that if the action of the 


1 C 1 -field is included along with the action of the 1 Av '-field, 
a universe will he obtained that approximates to the steady- 
state universe on a large scale although there may be local 
irregularities. In this universe, the value of 0 will be 
finite and its gradient time-like and of unit magnitude. 

Given this universe, we may check it for consistenc 7 f by 
ciaculating the advanced and retarded 'O'-fields and finding 
if their sum is finite. We shall not do this directly but 
will show that the advanced field is infinite while the retar¬ 
ded field is finite. 

Consider a region in space-time bounded by a three- 
dimensional space-like hypersurface D ^^ "tj he ^ ~ i 

and the oast light cone of some point r to the future 


of 




By Gauss's theorem 


Q C /- Q dx 



Zt-V 


2 CcL$ 

d a 



J K ; k /- 9 X 


Let the advanced field 

✓} t 1 Z'" 1 h 

Then L and 


1 / 


i P 1 
be U 


d£ 

d 


produced by sources within 

will be zero on , and hence 


re 
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X 




r 


J P 


Kcis 

a ix 


r K 

But is the rate of creation of raatter= hu (const.) in 

JK 

the steady-state universe, and hence 


mV 


D 


ru 


V. 


as the point is taken further into the future, the volume 

of the region 1/ tends to infinity. Rov/ever, the area of the 

hvpersurface 5 tends to a finite limit owing to horizon 
effects, 'therefore the gradient 


c) 


must be infinite 




A similar calculation shows the gradient of the retarded 
field to be finite. Their sums cannot therefore give the 


theory 


field of unit gradient required by the Hoyle-Narlikar 
Xt is worth noting that this result was obtained 
without assumptions of a smooth distribution of matter or 


of conformal!flatness. 
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5* Conclusion 


It is one of the weaknesses of the Einstein theory of 
relativity that although it furnishes field equations it does 
not provide boundary conditions for them. Thus it does not 
give a unique model for the universe but allows a whole series 
of models. Gle arly a theory that provided boundary conditions 
and thus restricted the possible solutions would be very 
attractive. The Hoyle-Marlikar theory does Just that(the 


requirement that dd 


) J_ y. 

(Yi . 't i rv\ 

let 


-A 




O 


equivalent to a boundary condition). Unfortunately, as we 
have seen above, this condition excludes those models that 
seem uo correspond to the actual universe, namely the 
Robertson-Walksr models. 

The calculations given above have considered the universe 
as being filled with a uniform distribution of matter. This 
is legitimate if we are able to make the ’smooth-fluid* 
approximation to obtain the Einstein equations. Alternatively 
if this approximation is invalid, it cannot he said that the 
tneory yields the Einstein equations. 

It might possibly be that local irregularities could make 
,VY ^ finite, but this has certainly not been demonstrated 
and seems unlikely in view of the fact that, in the Hoyle- 
warlikar direct-particle interaction theory of their 1 C-field, 
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which, is derived from a very similar action-principle, it can 
"be shown without assuming a smooth distribution that the 
advanced 'C* field will be infinite in an expanding universe 


with ere ation. 


file reason that it is poss 


to formulate a direct- 


particle interaction theory of electrodynamics that does not 
encounter this difficulty of having the advanced solution 
infinite is that in electrodynamics there are equal numbers 
of sources of positive and negative sign, i'heir fields can 
cancel each other out and the total field can be zero apart 
from local irregularities. This suggest that a possiole way 
to save the Hoyle-Narlikar theory would be to allow masses of 
both oositive and negative sign. The action would be 


4 


■r 

(a 


are gravitational charges aralogou 


s to 


where Q a, > c j , b 

electric char ges. Particles of positive <J. in a positive 
Vo ’-field and particles of negative <j_ in a negative 1 m 
field would have the normal gravitational properties, ^ 
tney would have positive gravitational and inertial masses. 
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A particle of negative 


in a positive ’/n '-field would 


still follow a geodesic. Therefore it would be attracted by 
a particle of positive O . Its own gravitational effect 
however would be to repel all other particles. Thus it would 
have the properties of the negative mass described by Bondi 
that is, negative gravitational mass and negative inertial 


mass. 


Bince there does not seem to be any matter having 
these properties in our region of space ( where M ^ const. > O ) 
there must clearly be separation on a very large scale. It 


would not be possible to identify particles of negative o, 
with antimatter, since it is known that antimatter has positive 
inertial mass. Hwever, the introduction of negative masses 


would probably raise more difficulties than itvould solve. 
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CHAPTER 2 


PERTURBATIONS 

1 • Introductio n. 

Perturbations of a spatially isotropic and homogeneous expanding 

universe have been investigated in a Newtonian approximation by 

f1 ) (9 \ 

Bonnor v and relativistically by Lifshitz v Liftshitz and 

Khalatnikov^ and IrvineTheir method was to consider small 
variations of the metric tensor. This has the disadvantage that the 
metric tensor is not a physically significant quantity, since one 
cannot directly measure it, but only its second derivatives. It is 
thus not obvious what the physical interpretation of a given 
perturbation of the metric is* Indeed it need have no physical 
significance at all, but merely correspond to a coordinate trans¬ 
formation. Instead it seems preferable to deal in terms of 
perturbations of the physically significant quantity, the curvature. 


Notation 


Space-time is represented as a four-dimensional Riemannian space 
with metric tensor g , of signature +2. Covariant differentiation 
in this space is indicated by a semi-colon. Square brackets around 
indices indicate antisymmetrisation and round brackets symmetrisation 
The conventions for the Riemann and Ricci tensors are:- 




; 0*0 


-r Oi 

— *L 


TV' 


o*c b ^ 




is the alternating tensor. 

Units are such that k the gravitational constant and c, the speed of 


light are one. 
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We will assume 


3• The Field Equations 

We assume the Einstein equations: 

± S.tR - - T„b 

where T > is the energy momentum tensor of matter, 
that the matter consists of a perfect fluid. Then, 

Tab - fA. U a U fe fi 

where U is the velocity of the fluid, U U a = - 1 : 

a ' a 


jj is the density , 
jx is the pressure 

* s the P ro 3 ec ' tion operator 
into the hyperplane orthogonal to U : 

EL 

h 0 t U 1, = o. 


We decompose the gradient of the velocity vector U_ as 

gL 


where U fl * 

B = 
<r ofc - 

b ^ 

flow lines U . 

a 


^ e< * b LI 


u*- 


U.(cjci) hot h b 

u t c i a 3 Hjl 


' 5 




is the acceleration, 
is the expansion, 

is the shear, 


is the rotation of the 


We define the rotation vector u)_ as 


oo 


<x 


Vua co tJ u b 


We may decompose the Riemann tensor R 


abed 


into the Ricci tensor 


R 


and the ’Veyl tensor C . ^ : 

R «i=c^ - Co.bc d - Rcj b * 3 Ri] 

Cabccl “ b] 


- R 




/-** Ok 

w \ U a. “ c 


C <x^ bed 


© Stephen Hawking, All Rights Reserved 












r is that -part of the curvature that is not determined locally by 

abed 

the matter. It may thus be taken as representing the free gravit¬ 
ational field (Jordan, Ehlers and Kundt^). We may decompose it 
into its "electric" and "magnetic'* components, 

E» b = - C. bp , , 

Hcb ~ >J<)rbs U p U ? 

c <* * i i t ^j £ Cp- 

tc*b * S u - A 8{v E ? 

-2 U p H nLt u d] , £^? cdrs u v H sL « u t] ? 

R a * - o 

E 0 ^ b = H= ° > 

E , and H . each have five independent components, 
ab ab 

We regard the Bianchi identities, 


R 


CkVLc.<ijS 3 


te o 


as field equations for the free gravitational field 


Then 


Cabcd' - c[b;aj + 'g' 



and Trttmper, °') 
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Using the decompositions given above, we may write these in a form 
analogous to the Maxwell equations. 


K* 

' E 


* C.<J 

h 

► 3 H a 

b 

.bto - 

^7 abeti ^ OT e 

H de = 

= y h a b 


(1) 

K 

b H 

be; ^ 

h * 

„ * P 

b ' 

y j 6 b c d ^ (T C £ 

i-de 

E - 

(/“ M 

to Q 

^ J 

(2) 

X 


ho' 

7 b) cd e 

u c H | 

E ab © - 

E (a 

c 



-- 


— 

E% 

- 

V occie 

^7 b p^f r U aP 

cT dt > 

E er 



&■ 

■ #’ 


+ A H d 

1 

« ^ bcefce 

u c u e 

* ^ t(^4) 

<Tab 

7 


(3) 


L cde u c H*v 6 “ H%W b)t 

, H (a. ^b) c“^a£<sle ¥ Jbf><jr'M'U* S 0* t! ^H / ^ 

f 2 H lJ a ^ be. ol e U fc U e C O 

where _L indicates projection by h^ orthogonal to U ft . 

{c,f, Trhmper, , 

The contracted Bianchi identities give, 

b ~ o , 

jx * ( w + fv) £ ^ ° > ( 5 ) 

(j* fx) U Q + J p>;b ^ ^ ” O - (6) 


The definition of the Riemann tensor is, 

a^lbe.j ~ 2 "RapU ^ 

Using the decompositions as above we may obtain what may be regarded 
as "equations of motion",- 
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Q = z Cjo 2- - 2 cr x - -f 9" + K 0 ;ft 


w , 


X GO a u, = - J ^ °c[ a ^j,i C + U C.P' si ^ ft ^ 


X <5~ 


a. b 


'.iJr 


ab ' CO^^CO b 


<3* C b ~ T^cb 0 


- -f h ob ( a to 2 - 2 .cr 2 + u c ;c ) + u a a b 


+ U(p^) hS hi 


where 2 co 2 - go a t cj aia , z <T - <-tob CT ab 


'Te also obtain what may be regarded as equations of constraint 


@;b H 


b 

a 


2 . 


( CO & c * b) h 


■*v 


U 


V? + 


( 10 ) 


<X> a 


r=. Z GO a U ^ 


H 


CM 


- h% ^b)ttie (__<■*)^ d;e t <7"f d '* j 


( 12 ) 


We consider perturbations of a universe that in the undisturbed state 
in conformally flat, that is 


C v , = 0 

abed 


By equations (l) - ( 3 ), this implies, 

O*o.b ” 00 ab - 0 

h«, b ^. b = O = 0; b ^ 
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If we assume an equation of state of the form, fv - j a (p) , 
then by (6), (10), f L; bh b « ^ o = ^ & 


This implies that the universe is spatially homogeneous and isotropic 
since there is no direction defined in the 3-space orthogonal to U , 

d. 

In this universe we consider small perturbations of the motion 
of the fluid and of the Weyl tensor. We neglect products of small 
quantities and perform derivatives with respect to the undisturbed 
metric. Since all the quantities we are interested in 'with the 
exception of the scalars, p, jx f 6 have unperturbed value zero, we 
avoid perturbations that merely represent coordinate transformation 
and have no physical significance. 

To the first order the equations (l ) - (4) a n( i (7) - (9) are 


E - r © b > 

H.t' k = u)» , 

# / {jjf ^ 

b ^ ^ ¥ IF Q ^ 6? t H £o fc>) c d e ^ ^ f — z 1 y^- 4 K) ^ at y 

H a b ¥ ^ab 0 ■* U<L ^ f. ' ~ o 7 

0 ^ - 5 " 6 4 + U-*' A - 9 

cj«!> = 6 * ^Lt:fl h « , 


cr.b = 






(13) 

(14) 

(15) 

( 16 ) 
(17) 


( 18 ) 


(19) 
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prom these we see- that perturbations of rotation or of or H gb do 

not produce perturbations of the expansion or the density. Nor do 
perturbations of E ^ and H ^ produce rotational perturbations. 


4. The Undistur b ed Metr ic 

Since in the unperturbed state the rotation and acceleration 
are zero, U must be hypersurface orthogonal. 

3. 


U 


*r\ 


T 

where 't measures the proper time along the world lines. As the 
surfaces T = constant are homogeneous and isotropic they must be 
3-Eurfaces of constant curvature. Therefore the metric can be 


written, 


d s 1 = - dr 1 + ,fT d y 


where 


a * ciw , 


T 


is the line element of a space of 


aero or unit positive or negative curvature, 


Vs define t by, 


then 


d t _ I 

Tr-n 


at S 


dv'3 


In this metric 


u 




“ (-O, o ; o,o) . 


Q ■ 


s£i . 5 a' 

a TP 

(prime denotes differentiation with respect to t) 
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Then, by ( 5 ), (7) 


P = 




‘7 


( 20 ) 


**■ 


3H 

a 




( 21 ) 


If we know the relation between p. and f\ , we may determine fl 
We will consider the two extreme cases, fa = 0 (dust) and 
(radiation). Any physical situation should lie between these, 

P or , j\ ^ 0. 

By (20), 


M 


= const 


•ri 

PL 




> n. _ 1 - 

M n 2 O 


- O 


_L — = f 

hA n ^ 


E 


const 


(a) Por E 0 


n = 


2F 


? ( 


Co 5 b tl t 


v 3 


(b) For E = 0, 


n - 




2 E 


(/X 

We/v\ 


iiVI 


t - 


t 


T = hi t 7 ■ 

c a6 L / 


/ 


(c) For E 0 


9 


Q ~ i_ £ ( { 


CDS 


EM 


t 


2 F 


t'-fe 


M 


sU V' 


EM 


E represents the energy (kinetic + potential) per unit mass. 

If it is non-negative the universe will expand indefinitely, other¬ 


wise it will eventually contract again 
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By the Gauss Codazzi equations *R,the curvature of the 


hypersurface T = const., is 


z E M 


a 


If E > o > 




* i 

■< * :ci a 

;- = o 
- ; 

K 

h - o 

E <° > 

r /' 

R " n? 

I 1 or ix ~ C /-5 

rrz-s-.- 9HSWIU ' ' AAMKmHHBII 




* 

i i 

-.ti 


/ 



(A - 


/v\ = 


* v 


? a 


(a) For E > 0 


j 


* 1 


* i 


3 Cl 
Ta 


pr <*i n h "£ 

£ 7 


X * 


^ cr-S-b t - | ^ 


a 


(b) For E = 0 


ft - t 


L. 


r\ - o 


(c) For E 0 


a = 


F 


SX H 


t 


J 


X 


CO'S 




Rot ational Pertur b ations 


-■* , UM 


By (6) 




k c , I- 

e*. ' 


1 b * 


CG e .ib ft 
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For “ft ® ° 


co •= 





For 






Thus rotation dies away as the universe expands. This is in fact a 

statement of the conservation of angular momentum in an expanding 
universe. 


6. Perturbations of Densi ty 

For ji-= 0 we have the equations, 

j, i ~ - ^.Q 

6 - -f e'-ir/A 

These involve no spatial derivatives. Thus the behaviour of one 

region is unaffected by the behaviour of another. Perturbations 

\ 

will consist in some regions having slightly higher or lower values 
of E than the average. If the universe as wh 'le has a value of E 
greater than zero, a small perturbation will till have E greater 
than gero and will continue to expand. It will not contract to 
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form a galaxy* If the universe has a value of E less than zero, a 
small perturbation can contract. However it will only begin 
contracting at a time 8 T earlier than the whole universe begins 
c ontrac t ing, where 



T 0 is the time at which the whole universe begins contracting. 
There is only any real instability when E = 0. This case is of 
measure zero relative to all the possible values E can have. 

However* this cannot really be used as an arguement to dismiss it 
as there might be some reason why the universe should have E = 0, 

For a region with energy -8E , in a universe with E = 0 






For E = 0, ju_ a ± 

/ -a 


r 3 

Thus the perturbation grows only as , This is not fast enough 


to produce galaxies from statistical fluctuations even if these 
could occur. However, since an evolutionary universe has a particle 



in the early stages. This makes it even more difficult for 


statistical fluctuations to occur over a rag:' y n until light had time 


to cross the region 
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r\< L * * 

•5 t/ *” jU ■+" W- a 


a. 


^ a “ 


- - h 


V, 




jb-i. 


as before, a perturbation cannot contract unless it has a negative 
value of Eg The action of the pressure forces make it still more 
difficult for it to contract. Eliminating 6, 


ay 


^ f 


* 1 






LA 


; Cl 


Cl 


- O 




^ iA ^ ]*} 


a b 


+" U a IA 


C4. 


J 

4* 






to our approximation 


h W tt Vb is the Laplacian in the hypersurface T = constant. 


he represent the perturbation as a sum of eigenfunctions S 


i(n) 


of 


this operator, where, 


s 


(*V) 


u c. o 


h"(Wt S'f 0 ; 


n 

Q 


s 


£rO 


These eigenfunctions will be hyperspherical and pseudohyperspherical 
harmonics in cases (c) and (a) respectively and plane waves in case 
(b). In case (c) n will take only discrete * .lues but in (a) and 
(b) it will take all positive values. 
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where ji is the undisturbed density. 



As long as /^ c > 




will grow. 




C x +- D V~ l 


These perturbations grow for as long as light has not had time to 
travel a significant distance compared to the scale of the perturbation 
( ^ Q ). Until that time pressure forces cannot act to even out 
perturbations. 

B % \ b'"° OL +■ ^ B w = o 

a 3 


When 


Y 1 


n 


» /A o 



’Ve obtain sound waves whose amplitude decreases with time. These 

(3) 

results confirm those obtained by Lifshitz and Khalatnikov . 

Prom the forgoing we see that galaxies cannot form as the resul 
of the growth of small perturbations. We may aspect that other non- 
gravitational force;.will have an effect smaller than pressure equal 
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to one third of the density and so will not cause relative perturbations 
to grow faster than T , To account for galaxies in an evolutionary 
universe we must assume there were finite, nan-statistical, initial 
inhomogeneities. 


7* The Steady-state Unive rse 

To obtain the steady-state universe we must add extra terms to 
the energy-moment inn tensor. Koyle and Narlikar use, 


^ M 


Ufe 


^ a.b O'«. 0 + /g. 3 o ^ 0 ©I 0 


( 20 ) 


where, 




y 

ii a. 




Ol 


0 


- (/*■ * t) w-1 

Since T^b^-o 

i- £/•+ ft) 6 -h u a C CL Cj, /lj o 



(21 ) 


( 22 ) 


There is a difficulty here, if we require that the 1? Cfield 
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should not produce acceleration or, in other words, that the matter 
created should have the same velocity as the matter already in 
existence,We must theft have 



However since C is a scalar, this implies that the rotation of the 
medium is zero. On the other hand if (23) does not hold, the equations 
are indeterminate (c.f. Raychaudhuri and Bannerjee- ; ). In order to 
have a determinate set of equations we will adopt (23) hut drop the 


requirement that C r is the gradient of a scalar. The condition (23) 


is not very satisfactory hut it is difficult to think of one more 


satisfactory. Hoyle and Narlikar 


( 12 ) 


seek to avoid this difficulty 


hy taking a particle rather than a fluid picture. However this has a 


o 


erious drawback since it leads to infinite fields (Hawking 


(13) 


). 


From (17), 

r~ 




AJ 



- (A 4 /0 - K) 0 


- 0 
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fl 


4" 


4 i 


f*- 4 /'v + tyQ J 


A j-1 j-C)Q 
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For 


(A 4 fc) =• 9 £ 1 ~ A 4 ftf) 


• (M- * fc ) 


i 


Thus, small perturbations of density die away 


Moreover equation (18) 


still holds, and therefore rotational perturbations also die away 


Squation (1 9) now beeomes 

0 = ~ 


3 


+1 




9 


4 >( 


to 


These results confirm those obtained by Hoyle and Narlikar 


(i 4) 


We 


see therefore that galaxies cannot be formed in the steady~state 
universe by the growth of small perturbations. However this does not 
exclude the possibility that there might by a self-perpetuating 
system of finite perturbations which could produce galaxies. 


(Sciama 


0 -0 , Roxburgh and Saffman ! ° ) 


8. Gravitational Waves 


rTTWm II j 


-t.v-JL-lAsZTlk-j r r '■ i”.~ 


We now consider perturbations of the Weyl tensor that do not 
arise from rotational or density perturbations, that is, 


E 


C J 


H 


a 


h 


“ O 


Multiplying (15) by U l V< 


fits 


b^ 


and (1 6) by h £ > ) ^ u r 
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we obtain, after a lot of reduction, 



In empty space with a non-expanding congruence U a this reduces to 
the usual form of the linearised theory, 

□ 4 - o 


The second term in (2h) is the Laplacian in the hypersurface 


"Tr = constant 


, acting on . We will write as a sum of 


ab 


eigensfunctions of this operator. 


£ 


o.b 


= Z 


ab 


where 


W 00 

V oV> = ° > 


(Veot;, hi k S b = 


v ^ 


\/ ^ 
v D i> - o 


V 


ot- 


- O 


Then 


E 


a. 


T A /Crt> x/ ^ 
2_ 


E 


a 




z 


A 


aV") — f 


a. 


t V _n 

n'‘ 


v 


O'!' 


Ot 


© Stephen Hawking, All Rights Reserved 

















Inches , 1 

LUlxii 


T 


cm 
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Similarly, 


6 " 


Oft') , , (n) 


C*. 


Z D' 'v ob 


Then by (19) 


fi-i) 


C2 A 


iri 


2 V 


Q 


Q 


Substituting in (24) 




6 Cl' 


A' M f A l " 1 


* -.n* ( o! x 

n + 3 4— r © 4^-r 

a n 3 - 


* D 


c*> 


r n’(j**h) ixQ - 0 


+ 3 


We may differentiate again and substitute for D* 


For vi»t 


and 


a » 






e 


t«t 


so the gravitational field E b decreases as Q" 1 and the '''energy” 


i(S ab E ab + H ab H ab ) as ft 


■6 


. We might expect this as the 


Bianchi identities may be written, to the linear approximation 


n c*b C oO = J.b6. 

o 


Therefore if the interaction with the matter could be neglected 


C . , would be proportional to O. and E . , H , to 

aoca ab ab 


•1 


In the steady**state universe when p and 0 have reached their 


equilibrium values. 


a 




* + h 


V 



\ J 


Cu 


be 
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Thus the interaction of the "C" field with gravitational radiation i 


1 s 


equal and opposite to that of the matter. There is then no net 
interaction, and E &b and H decrease as Q. . 


The "energy" -|(E, 


,-ab „ v „ab 


ab 


+ H. b ir*“) depends on se*ond derivatives 


of the metric. It is therefore proportional to the frequency squared 


times 


the energy as measured by the energy momenturn pseudo-tensor, in 
a local co-moving Cartesian coordinate system which depends only on 
first derivatives. Since the frequency will he inversely proportional 
to O , the energy measured by the pseudo-tensor will be proportional 
to CL ^ as for other rest mass zero fields. 


9. Ab sorption of G ravitational Waves 


As we have seen, gravitational waves are not absorbed by a 
perfect fluid. Suppose however there is a small amount of viscosity 
We may represent this by the addition of a term "X to the 

energy-momentum tensor, where "X is the coefficient of viscosity 
(Ehlers/ 17 '). 


Since 


— o 


we have 


f*- * {ft* tO 9 


> <T L - 


(25) 



fib ^ 


X * > <n v ;b ’ o 


( 26 ) 
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Squations ( 15 ) ( 16 ) become 

+ hV?b)cole ut " € - ' -k fa + fv.) <J^b 

-4 X ( E«b / (27) 

■'Lb 4 H oi) 0 ~ H V (ec7b)ccU U C E^" e - -J, X Heib ( 28 > 

e^-B, ■ 

The extra terms on the right of equations (27), (28) are similar to 
conduction terms in Maxwell's equations and will cause the wave to 
decrease by a factor £ . Neglecting expansion for the moment, 

suppose we have a wave of the form, 

P . * F P ivT 

Thi 3 will be absorbed in a characteristic time 7 a independent of 
frequency. By (25) the rate of gain of rest mass energy of the 
matter will be which by (19) will be 2.X Jo V . Thus the 

available energy in the wave is 4 1 ^ • This confirms that the 

density of available energy of gravitational radiation will decrease 
as Q in an expanding universe. Prom this we see that 
gravitational radiation behaves in much the same way as other 
radiation fields. In the early stages of an evolutionary universe 
when the temperature was very high we might expect an equilibrium to 
be set up between black-body electromagnetic j diation and black-body 
gravitational radition. Since they both have two polarisations their 


f 
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energy densities should toe eaual. As the universe expanded they woul 
both cool adiabatically at the same rate. As we know the temperature 
of black-body extragalactic electromagnetic radiation is less than 
5°S , the temperature of the black-body gravitational radiation must 
be also less than this which would be absolutely undetectable. ITow 
the energy cf gravitational radiation does not contribute to the 
ordinary energy momentum tensor T . , Nevertheless it will have an 

active gravitational effect. By the expansion equation, 

6 . 


for incoherent gravitational radiation at frequency v , 



r-*. 

E v 



But the energy density of the radiation is 


t- 2 - ' z 

4 E v 


V * 


e 


5 @ ' 'Z/ u -Q - •" 3 fa ) 


where is the gravitational ’’’energy'* density. Thus gravitational 

■■ ■! 

radiation has an active attractive gravitational effect. It is 
interesting that this seems to be oust half that of electromagnetic 
radiation. 

It has been suggested by Hogarth! 1 ^ and Hoyle and Narlikar^ 1 , 
that there may be a connection between the absorption of radiation 
and the Arrow of Time, Thus in universes like the steady-state, in 
which all electromagnetic radiation emitted is jventually absorbed by 
other matter, the Absorber theory would predic retarded solutions of 
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1,1 1, 

1 i i 
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> 1 1 
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cm 
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3 

I | i 

4 1 \\ 


the Maxwell equations while in evolutionary universes in which 

electromagnetic radiation is not completely absorbed it would predict 

advanced solutions* Similarly, if one accepted this theory, one would 

expect retarded solutions of the Einstein equations if and only if all 

gravitational radiation emitted is eventually absorbed by other matter* 

Clearly this is so for the steady-state universe since X will be 

constant. In evolutionary universes X will be a function of time. 

We will obtain complete absorption if diverges. How f <■ ? a gas, 

i ^ 

NocT 2 where T is the temperature. For a monatomic gas, T , 

therefore the integral will diverge (just). However the expression 
used for viscosity assumed that the mean free path of the atoms was 
small compared to the scale of the disturbance* Since the mean free 
path oc [i oe il and the wavelength Cl “ , the mean free path will 

eventually be greater than the ’wavelength and so the effective viscosity 
will decrease more rapidly than Q _1 * Thus there will not be complete 
absorption and the theory would not predict retarded solutions. 

However this is slightly academic since gravitational radiation has not 
yet been detected, let alone investigated to see whether it corresponds 
to a retarded or advanced solution. 
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CHAPTER 5 


Gravitational Radiation In An 


expanding Universe 


Gravitational radiation in empty asymptotically flao 


space has been examined by means of asymptotic expansions 


by a number of authors. 


( 1 - 4 ) 


•They find that the different 


components of the outgoing radiation field "peel off", that 
is, they go as different powers of the affine radial distance 


ea 


Ii one wishes to investigate how this behaviour is modifi 
by the presence of matter, one is faced with a difficulty 
that does not arise in the case of, say, electromagnetic 
radiation in matter. For this one can consider the radiation 
travelling through an infinite uniform medium that is static 
apart from the disturbance created by the radiation. In the 
case of gravitational radiation this is not possible. For, 
if the medium were initially static, its own self gravitation 


would cause it to contract in on itself and it would cease t 


o 


be static. 


Hence one is forced to investigate gravitation 


O-- _«_ 


radiation in matter that is either contracting or expanding. 
i j .s in Chapter 2, we identify the Weyl or conformal 


tensor l v .bcc L 


with the free gravitational field and the 


9 


Kicci-uensor K,^ with the contribution of the matter to the 
curvature. Instead of considering gravitational radiation in 


f 
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space that approaches 


asymptotically flat space, tnat is, 
flat space at large radial distances, we consider it in 
asymptotically conformally flat space. As it is only 
conformally flat, the Ricci-tensor and the density of matter 
need not be zero. 


To avoid 


essentially non-gravitation&l phenomena such 


as sound waves, we will consider 
travelling through dust. It was 
conformally flat universe filled 


gravitationalradiation 
shown in Chapter 2 that a 
with dust must have one of 


tile metrics 
(a) 


cl s 2 . (ct t 1 - cl f 

SI , R ( 1 - c ) 








Type(a) represents a universe in which the matter 
expands from the initial singularity with insufficient 
to reach infinity and so falls hack again to another 


P in 

^ -*■ a>d 


singularity. 


It is therefore unsuitable for 


& 


H "3 

-i- 


seussion ox 
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since one cannot 0 et an infinite distanc 


technique 


gravitational radiation by a method of asymptotic expansions 

oe from thi^- source. 

i./P® Ch) is tne binstein—i)e ditter universe in. which the 

mailer has just sufficient energy, to reach infinity. It is 

ohus a special case. f. Norman has investigated the 

"peeling off" behaviour in this case using Penrose’s conformal 

( 6 ) 

* He was However forced so make certain assumpt¬ 
ions about the movement of the matter which will be shown to 
be false. Moreover, he was misled by the special nature of 
the Hinstein-De Sitter universe in which affine ana luminosity 

f y 

distances differ. Another reason for not considering radiation 
in the dinstein-ne Sitter universe is that it is unstable. 

lie passage of a gravitational wave will cause it to contract 
again eventually and develop a singularity. 

.0 v/ill tnerefore consider racxiation in a universe of 
type (c) which corresponds to the general 


case wnere the 

matter is expanding with more than enough energy to avoid 
contract!ng again. 


ihe Newman-Penrose Formalism 


we employ the notation of Newman and Penrose,^ A 
tetrad of null vectors jiT fl 1 ^ 


n 


is introduced 



















o 



we label these vectors with a tetrad index 


'Z 

(X 




= /V 




; 1 




O 

5 


<* . 


tetrad indices are raised and lowered with ohe metric 


Qibi 


rie have 


ib 

n - 
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G O 
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1 o o 


0 O O " j 

0 0—1 0 


/-i _, v 

2 2 
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M . v/ ^ u- , v . M ~r~ ^ 

i n + fi l - w fl) - 


— M l' 

/r /n 
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Ricci rotation coefficients are defined by: 


be 

V 


\\ 


0 1 b 
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0 - 
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4 t 
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2 ^ 


*/ 


( 2 . 3 ) 
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In fact it is more convenient to work In terms of twelve 
comnlex combinations of rotation coefficients defined as 


follows: 
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3. Goordinates 


Like Newman and Penrose, we introduce a null coordinate 
oL ( = X ' ) 

^ Vfi “ «o (LO 


A 

nj 


tX , i, LL ■ v 

) jn 1 


we take 


a 


• Thus L a 
geodesic and irrotational. ’Ibis implies 

T - o 

e ■- f 

i , -€ 

cL " r ^ 


will be 


*■) 


we t ake 


fit AvG n 


L 


h 


to be parallelly transported 


along [_ . This gives 


I 1 


£ 






< 

As a second coordinate we take an affine paramete 
long the geodesics L aa 

K L r * / 


r/-. xb 


(Li) 


j i •/* 

X and X' are two coordinates that label the geodesic 


in the surface Li = const 

x J / ^ • 

iui/ A ;/ u L ■ 


Thus 
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i I 


\ 0 l DO 

I i f 1 f ‘ T’ 
o f ! S y T < ' 
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In these coordinates 
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The Field Ecuations 


We may calculate the Ricci and Weyl tense 


r components from 


the relations 
c,cL ^ i? c; fl. 
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Expressing the rotation coefficients in terms of the metric, 
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As in Chapter 2 we use the Bianchi identities as 
field equations for the Weyl tensor. In the hewman -Penrose 


formalism they may be written: 


(I am indebted to R. G. McLena^han for these) 
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4, The Undisturbed Metric 


The undisturbed metric may be written 

ds* Stride 2 -- cl (5 * - s i/v/dp C td 2 >- 5 i^d cp^ 


put 


hen 


SL - 0 C cos A t - 0 

Lt - *> -p 

• ’ 

_ SL 1 r-ctiA Z -r ScLuoit - S//V 



ly_ is a null coordinate 

To calculate , the affine parameter, we note that C- 


is an affine parameter for the metric within the square 

- £<v,& f) 


brackets. Therefore 


will be an affine parameter fori'*/) 


cU- 


(jL± ) 



is constant along the null geodesic, No finally it 


would be taken so that Ned) when c ; 1 


However, 


in our case it will be more convenient to make it zero and 


define 


r 


as 


r = pJl^oO 

J o 


£ <i- 3) 


This means that surfaces of constant (r are surfaces of 
constant £ . This may seem rather odd, but it should 

be pointed out that the choice of £ will not affect the 


asymptotic dependence of quantities. 

f , o(r~ rt ) 

Than v 

^ oir " fV ) 


That is, if 


f 


% 

f 


: ffS 
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It proves easier to perform trie calculations with this 
choice of h but all results could be transformed back 


to a more normal coordinate system, 

(U-V ^n 


Erom 


f _ ff x /-Lsi/vk^ ~9St*~k hr r 


(b. t /) 


The matter in the universe is a. 




ssumed to be dust so its 


energy tensor may be written 

Tst = h ^ V b 

i?or the undisturbed case, from Chapter 


4 * -> 3 




V 


£ R 

IP 

Sic. 


o- 


1 / V 0 -* I 


/ U 6 ) 


N o w 


JfL 


S 79 


lof\$ 
S 


o (S'") 


&S) 


wnere 


S’ 


1 


r 


Therefore if we try to expand 
the result will be very messy and will 


as a series in power of I" 


involve terms of 


the form 


s' 


1 CL 


*It should be pointed out that the expansions used will 
only be assumed to be valid asymptotically. They will not 
be assumed to converge at finite distances nor will the 
' Quantities concerned be assumed analytic, (see A. Erdelyi 


Asymptotic Expansions - Dover 


> ■ 


© Stephen Hawking, All Rights Reserved 



















This does not invalidate it as an asymptotic expansion hut 


it makes it tedious to handle. For convenience therefore, 
we will perform the expansions in terms of ^ r ) 


will be defined in general as the same function of f* 
it is in the undisturbed case. That is 

t'i) 

, fp- fi sudifc- t r %_t 


whi c h 


as 


waere 


then 


1 


cl ■- 






d, r 




J1 


_ n 




/ 9 5 y 

£ 2 .* ^ 


ji 


For the third and fourth coordinates it is more convenient to 
use stereographic coordinates than spherical polars. 

Since the matter is dust its energy-momentum tensor 
and hence the Ricci-tensor have only four independent 
components. We will take these as 4,<L.4c 
(Since 0> O( is complex it represents two components) 


In terms of these the other components of the Ricci-tensor 


may be expressed as: 

id * 


% 


4 


O i 


4) 


90. 




oo 


4 > 




AM 

$0* 


i + 



-1 
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For the undisturbed universe with the coordinate system 


given: 













Using these values and the fact that in the undisturbed 


universe all the 



are zero, we may integrate equations . 


(3. 10 - 50 ) to find the values of the spin coefficients for the 



r " 
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5* Boundary Conditions 


We wish to consider radiation in a universe that 
asymptotically approaches the undisturbed universe given 


above 


<t>„ 


and 


A 


will then have the values 


.riven above olus terms of smaller order 


To determine thi 


s 


order and the order of 




and 


% 


, there■are two 


ways in which we may proceed, be may take the smallest orders 
that will permit radiation, that is * Qf \T y 

<L 


Larger order terms than these in 


, A 


ana 




turn out to have their 


CL 


Pi 

dependent 


derivatives 

-I 


of 




only on themselves and not on the r coefficient 

, the radiation field. They are thus disturbances 


not produced by the radiation field and will not be considered 


Alternatively we may proceed by a method of successive 
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a 


pproximations. We take the undisturbed values of the 


soin 


coefficients ano. use then to solve the Bianchi Identities a 


boundary condition that 

. , i/n f 

t nese 


■b 

Then 


the flat snace 


field equations for the conformal tensor usins: 

Tf a ' , oh '/ 1 

in equations ( 3JO - ) calculate the disturb- 


substituting 

w 


r 

ances induced in the spin coefficients and substituting these 
back in the Bianchi Identities, calculate the disturbances in 


the 


further iteration does not affect the orders 


of the disturbances. 


aoth these methods indicate that the boundary conditions 


should be: 

i\ 

o 

' 

o(si- 1 ) 

(hi) 

^OO 

- 3 R 4 . 

^ " J f 

SI* 

0 (SlJ 

fti) 

fo t 

~ 0 (SlJ 

(see next section) 

(i j) 

Jo 

-- 0 (-Jl- 

i 


(M) 


We also assume"uniform smoothness", that is: 


X 1..x * - ocsij 

J_A - - V) * ofJl'*) 

(,Sl u 

0 tj C # * 4 / 
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Inches 


cm 


2 


van be shown that if these boundary conditions 
hold on one hypersurface (u. = const.) they will hold on 
succeeding hypersurfaces and that these conditions are ■ 
most severe to permit radiation. 

Integration 

As Newman and Penrose, we begin by integrating the 
equations (5. 10 &11) 



where 





r 



Let 
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then 


pp P ■ 


let 


? 


j— CD 

-0'l) V 


then j) ^ y - — Cp Y 

since f >* (jpcLr o£2 



However 


Pv - 

V , 

<e - 

therefore J) 2 - V 
therefore J) v j 

V 

V 


P T Oft) 
rf t- OCr) 

o(r 'I) 

„ - ?- cp/- r O (y 

z (t 2 ) 


r T r- ofr *) r t t 

~r " ‘X *r o (r *) 


where /" is constant (&' 


1 


) -- c (r 


I 

a- 



S') 


is constant 


Li) 


ly** 

if (j- is non-singular (The case r singular corresponds to 
asymptotically plane or cylindrical surfaces and will not be 
considered here), 

"r o(r ) ~-2 f).~ 3 ' r c(sf) 

Q <5 “ (‘J) 

3 3 ). 

k S). '* 


Thus 


Let 


f 

G 

f 

6’ 


y 


O Cr'i) 

- 1 Jl' 1 


o(sr^) 

$ 


d-i) 


wnere 


- o6) 


Then using: 

D 1 


sr'( it 


v)y 


1 


R l n[ 

2 


i I 


hji 
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Integrating, 

3 


1 

Jli 0(0 



therefore 




For rL 
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IL (jlrOO)) ; ~Lr 0(SL'*) 

)Sl K > 

therefore 

L , Q(Sl_'‘) (jJQ\ 
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Repeat the process with 


p - 

- 2 fiSl'K 3 J1' V 

(b. lUj 

6 - 



where ^ - 



k. -- 

o 6) 


then c) ^ /j 

1-r 6(i) - 0 (si ' ') 


Cil 1 

K ; 6 B ( a ; x c ). r o(Sl ') 

pjj) 

<hi 

Sl-r - O (S)^ 1 6oc.SX^j 


kJ / 

- o v 0 C'^ 6o^Sl.^ 



Unlike Newman and Unti, we cannot make 


the transformation Y * & ~ p 


alter the boundary condition 


A 


£ 

liSl- 



zero by 


since this would 



Continuing the above process we derive: 

( „ -251 - Rsi'% f> Si V 
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cm 
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To determine the asymptotic behaviour of y., ■* £' s 


and oJ> 


nr* 


vie use the lemma proved by Newman and Penrose: 

matrixB&nd the column vector O are 


The f\ x f) 


given functions of x such that: 

- 2 


B ■= o(x 


b - °(x 


- 2 


(Ufi 


The n x n 


matrix fr is independent of x and has no 


eigenvalue with positive real part. Any eigenvalue with 
vanishing real part is regular. Then all solutions of: 


X ‘ T 2* ) 


are bounded 
vector. 


p) oD 


is a column 


>or reaso s to be explained below, we will assume for 


the moment that 


o) 


5 


= 

<fo, Z 





Is 


< I 


-S 


O (51 


-<» 


4 


O i 


/ 






We take as 




the column vector 


£il‘ •*, -n* », *, ^ fs.-n f W# 11 Hi 


/' o n J 

C. a ay 
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By equations 3. *5 / 


Pr - 


i.u, 

/ 


I'llUS 


. 0- 


and 

S'. 


b 





are 
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-i 
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<551 


0 1 


1 

"\ ‘ 
OK 





lO 


■ ? ) 

= oCSL*) 
= 06 ) 



- 3 o 6 R bft o o o o C i 

O O ° o O o O o o O 
0 
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Oqo aoooo o o 

OOH -I 0 o 0 o 0 “2 


3.43,3 


O I 
o 
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I 0 O -ioo o o 0 


0 i 

oj 

o| 

o 
o I 

I 

oj 

•2 


i, 7 t) 


expressions involvin ; 


3 d> 


<551 


i°‘ 



. 2 




^mce 


'■c 


= ol r!2 


X -- 
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Jsing this we integrate equation 


(3 . / i) 


by the same 


method as abov 


r .i 

w * 


We obtain 






X 



opv 3 ) 


r 


A. 23 


«'e may make a null rotation of the tetrad on each null 


geodesic 


u. 

a 7 yu 

A|> 


L jU 

rtyU 


a 


/- 


<X A"! -jj. r ^ ^ f^~ 


H 




r a 


•, £ 





is constant alone the geodesic since the tetrad is 


paralleily transported. 


3y taking 


i 


— > ^ f 
CJU - 2 L 


we may make 


* a 

a- - w 


Under a null rotation 

f 0 . : to, r “to 

Thus until we have specified the null rotation we cannot 


imoose a boundary condition on C) 


o f 


more severe than 


t 



We will specify the null rotation by 


a 


O 


and in that tetrad system will impose the 
boundary condition that <p Q - 0(SI /and is uniformly smooth. 


Then by using this condition on <J 


ana 


by equation 


(i-***) 


u 


■- o(dT') 
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using this in equation 


O 


then by equation 


^ O (SL 


- w 


L - 


putting this back in equation 

— O C 

by equation 

^ * 

by equation (?. 


(5. </ <f) 




o£Jl ^ 




by equation 


(j.ov) 


uj = LO 


ojT" 2 ?- 0 fJl) 


by equation 


o(^ 




by equation J'2-) 

by equation (svv) 


r 


^ _ o 


L 


o 


-2 


(J ^ 


r 


0 


(6 ■ 25~ ) 




6 - V 


By differentiating the equations used with respect to 

u- 


ono may snov 
uniformly smooth 
Adding equations 


w that r> y 


uO 


3 ^ 


'\l/ > 


i^u* ^ =■ ^ ^/o 


are 


and: 5 * & O ; - 
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'X 


i o 


T 


"t 
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we may use trie lemma again with 




1 


si* 

si 2 

[-h 


i 

T: 

X 


n 

x 


i 


fu~ 


By equations 3 # /£ 3 , /> ; £. 2.8* 


A 


1 - 2 - 0 


O 


£ and are 


0 


o 


o 


o 


- i 


Therefore 


A 



) 


J 


£ 


L ^ 


0 (J~L~ *<) 

ofJi'V 


) 





J 



and are uniformly smooth 


From f& 2 c i) and 


h 


using this in 


(2.1?) 

£Sl~ 

2 

(&. is) 


we may show 

-2 


l 


~T 


0(J\ " (a 


J 





0(SL % k>$Jl) 


then hy 


&• *) 


h 


£fi-‘r 




x L )Jl 


-2 


r 



^ he X]S) 


and 
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Integrating the radial equations 'X. / 3 X / 


1. C' £ X y S' 
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^si^^siKUzr 


V- f * ° r * °p 

2 -Qp°Jl- \i(z f)\ /O °lS 0 -r J2°6" °' 


r i (l °~ fi °/S 

■ 4 r * J ) 
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'-f* 

^ _ \° Pi '2. , / \0 dh\ 


A AT i-4 <A\ 
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- v 


.no 


lo /? nc> 

-4 6 * 




i f a - c> 

uj + r 
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if Mr J'V 
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-5 
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^4 

(£• ?j ) 


■rj 


O 
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^ j -XX 
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V 0^7 


p f 


)U 


-s 
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Adding equation 
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t 2x (i. s’9 ^ 
S 6 _ r 21 /) 


TlO 


2.>f r 2 
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— oC 
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'Therefore 


By equation 


ti = 

(2- it) 

i/ - 1 / 



<P 


o 



- V 


t - 


0(S1 


o 


X 


*\lr 


O 



-2 
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By the orthonormality relations 

^ y L it L ~ 

= A " \ cJ t 





5 


3 




3 


j- 






tor ,o 

a W Y 


ff 





‘J Yfl- 


By making the coordinate transformation 


5* 



- 1 


) 


4 < 


L> 


i 





'<L12> 



■ U°) 


i.u 



) (i-AP 


LL 




1 

f is 


U- 

r 

V 

t 

X T 


c L U, 


i 

V 


waere 


o 


e 

c 


X 


30 


I r 


c\)~ 






-x^( 







c 


J y 


' * C ’) - A ! “C 


X" = O 

We still have the coordinate freedom 







We may use this to reduce the leading term of <o ^ (^/ 
a conformally flat metric (c.f. Newman and Unti), that 

m _ 5?pyV/urh 2/ui'Xh h-rr9 


AO 


to 


4 I 


i 


is : 

Os 
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wnere 


30 


o 


L 


= PA, X U ) 


7 • Non-radial ‘Squations 

By comparing coefficients of the various powers of UL. 

AA 

in the non-radial equations of t 3 


, relations between 


the integration constants of the radial equations may be 


obtained 


In equation 


' 3 . 2.3 




the 

-°) 


rerm in 


si 


lR 


therefore 


M 






f 


- I 


therefore by 


(l ■ z y) 


u 


i n 1 


U f- 


0 (- 2 “‘) 


, %■ 

In equation AAS’oJ t the constant term is 


therefore 


\/° 

</°= 0 


is 


> 1 ) 


h_ 2 ) 


( 7. S 


Z 


By the jl tern 


m in {2. 


By the A term 


P~P. - 

J I 



C 


p - p. - 0 

f ¥ J t 


tiler of ore 


~ o 


r 


)? 


Jr 0 -- 0 


<z_Jd 
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if 


s 


iSie 




JJ 


y making a spatial rotation of the tetrad. 

- L <1 

, e m 


pfr .. 

we make s real. We take 2 : 


^2 





stereographic project factor for a 2 -sphere 


By the JTL ^ term in equation 


(IMS’) 


oL 


o 


-^V5e 


LA. 


where 


V 
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-r 




By the jV ' term an 

,ie* u ($VV Sr 

o n 2 


G-2-0 

sy vs) 
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ir 


jj, - -/r~ V V^Se 
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. V 

(S.60) 4.6 0 






[/ 
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where 


H 


ft 




th 


30 


O 


*1“ 
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ft 


2. ^ c 
e 




By making a 


? q 

coordinate transformation* 

t ' = f- - H dU^ 


S' 


// 


OlO 


ti ' - O 


therefore 


u 


, Bft-Qlsft- 


(d 


6 


€ 


By tae WT 4 term in A. 26) 

■ -fiftlit- e ) 

a 


_ ^° T LO 




therefore 


P„ 


?- O A 3. 


3 


-+* 


6 


M*- ±ft*t 


P 


o 



2 u. 

e. 1 *'■ 


Cf*‘ 


2 






*This transformation does not upset 

<7 


the boundary conditions 


on the hypersurface 


0^ — 
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?. 2 S') 
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By the jfi 
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By the ^Q/ 2 - term in 

to I 


0 - 22 ) 

(s l VS 


a £ 
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9 (h-€> 


(3. So) 


L^.2 
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. U 2 
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(l - 2o) 
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e u s 


therefore 


the refore 


O 


UJ 


- “ 2a' 


a T ■ 
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^ 1 -Ah 
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Using 


(f.cn.ih), 


t 


m 


ChJ2) 

;> . / ? ) 

C ^ • /£) 
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1 /?) 

J A Jl~'-A}-(i +t, 

<2 v 

+ ofji" 1 Jtj 
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2 t. 


F’ # 
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Using this in 
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0 fJL' 0 
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Cl. SI) 

3 
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3 

A 

A. A 


S' (). 60) (f. 6 /) 

= 0 CSV*) 

^ O (si'*) 

, o fil'") 


,) 


cv W 


A 

cl i/: 


A 


L N. 


O 


o C-^l -9 ) 


l-tvi) 


hold on 


Therefore if the boundary conditions ( A. / 
one null hypersurface, they will hold on succeeding hypersur¬ 


faces , 

By ($■ S?) ^ 

A - 0(01 2- 



"peeling off' 1 behaviour is therefore: 
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As mentioned before 
of che zero of r 


, this asymptotic behaviour is independent 
and will hold for any affine parameter r 


To perform the remaining integrations we 


will -assume 


for definiteness: 
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Thus che O 


derivative of dh ' , depends only on itsel 




and not on the radiation field. 


it therefore represents a 


type of disturbance unconnected with, radiation. If it is 
aero on one hypersurface, it will remain zero. In this case 
it is possible to continue the expansions of ail quantities 
in negative powers of S2. without any log terms appearing. 


6 . The Asymptotic frouo 


The metric has the form: 


tf 


i 3 


1 i 

• 2.1 2 . 


Jl 




- - o 

i ■ •* 


% t 


3 


'if 


o (si "*) 


3 


I'Z. 


I 

! 


-2P' X S ^ JZ'V ,. 2bp‘ z S o(2l 


— £? 


J 


The asymptotic group is the group of coordinate transformations 
that leave the form of the metric and of the boundary conditions 
unchanged. It can be derived most simply by considering the 
corresponding infinitesimal transformations: 

*._-J 
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ds 


J I 
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2 . 


cCfc k * i ^ /f 

' *■ / 3 3 ^ 
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3 lU 



To obtain the asymptotic group we demand 
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Inches 


r 

T 

T r T l i 


cm 

1 

\2 |3 

4 


reduced the leading terra of 4 to a conformally 


form. 


hue the only allowed transformations of x 


are t ne 


conformal transformations of the form: 


3 « w 

X "f LX 




C)sJL±L jj /? 


C (X 3 t lx* 


*r 


a 1 


Ac 


i 

/ 


then the six parameters a. ( A d, ci are given K* is uniquely 
determined hy (*<?”. /d) K^ is also uniquely determined. 


T hu a 


the asymptotic group is isomorphic to the conformal group in 

(?) 


two dimensions. 


has shown that this is isomornhic 


to the homogeneous Lorents group. 


It is also however isomorphic 


to the group of motions of a 3-space of constant negative 
curvature which is the group of the unperturbed fiobertson- 


v.alker space. Thus the asymptotic group is the same 


xs the 


group of the undisturbed space. It is not enlarged by the 
presence of radiation. This is interesting because in the 
case of gravitational radiation in empty, asymptotically flat 
space, it turns out that the asymptotic group contains net 
only the 10 dimensional inhomogeneous Lorentz group, the group 
of motions of flat space, but also infinite dimensional 


II 


supertranslations 1 '. It has been suggested that these 

supertranslations might have some physical significance in 
elementary particle physics. The above result would seem 
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Inches , 

lS i I i l 


7 

cm 1 1 


to indicate that this is probably not the case since our 
universe is almost certainly not asymptotically flat though 
it may be asymptotically Kobertson-walker. 


9 . VJhat an observ-er would measure 


The velocity vector V of an observer moving with the 

*Y\ 


dusL will be 


-S' 


jy' T o(ai~y 

ijl + 0 fJV 3 ) 


O (Sty 

0 (SI ' 3 


, the projection of the wave vector 


/ ° L ' 

C in t 


observer's rest-space ( the apparent direction of the wave) 


will be 


ha 


_ / v 




ho * 


~ a 
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T 


0 (SI 


- 4 



© Stephen Hawking, All Rights Reserved 



















o C ~ u ) 

0 (SL">) 



The observer’s orthonorraal tetrad may be comoleted by two 


soace-like unit vectors 


S 5 
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and C 
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- /C< /-J 

= ( ^ , <1, S 


(X 


o 


3-w measuring' the relative accelerations of neighbouring dust 


r J 


particles, the observer may determine the ’electric* 


components of the gravitational wave t 
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CHAPTER 4 


S in ru 1 ari t i e s 


are 


If the Einstein equations without cosmological constant 
.atisfied, a Robertson-Walker model can 'bounce 1 or avoid 


a singularity only if the pressure is less than minus one- 
third the density. This is clearly not a property possessed 
by normal matter though it might be possessed by a field cf 


negative energy density like the 'C field. However there is 
a grave quantum— mechanical difficulty associated witn o’ne 
existence of negative energy density, for there would be 
nothing to prevent the creation, in a given volume of space- 
time, of an infinite number of quanta of the negative energy 
field and a corresponding infinity of particles of positive 
energy. If we therefore exclude such fields, all dooertson- 
Walker models must be of the 'big—bang* type, -hat is taey 
have a singularity in the past ana maybe one in ijhe j.utu. 
as well. It has been suggested that the occurrence o 




these singularities is a consequence of the high degree of 
symmetry of the Robertson-Walker models which restricts the 
expansion and contraction so that they are purely 

■A? 

uhat inoi ?0 2 ?e al istic niod 0 1 s with f 0 w62? or no eXctC v s y mi.© 11 i- 
would not have a singularity. This chapter will be devoted 


and 
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to an examination of this question and it will be shown that 
provided certain physically reasonable conditions hold* any 


mo 


del must have & singularity, that is, it cannot be a 


^ 2 

geodesically complete C , piecewise G manifold 


2. The iundamental Equation 




The expansion 


0 * 1 /. 


; O- 


congruence with unit tangent vector 


of a time-like geodesic 
1 / obeys equation (7) 


of Chapter 2: 


B 1 / 




- ±B Z - L t % 

a u 


t 


- U v ' v 


CO 


A point 




A will be said to be a singular point on a geodesic 
of a time-like geodesic congruence iftffor the congruence 


is infinite on 


. A point Q will be said to be 

/ if it is a 


Y at % 

conjugate to a point P along a geodesic j 
singular point on V 1 of the congruence of all time-like 


geodesics through 


p . A point 0^ will be said to conjugate 
to a space-like hypersurfs.ee. b if it is a singular point; 

of the congruence of geodesic normals to fy .An alternative 
description of conjugate points may be given as follows: 


let 


Ov 


be a vector connecting points corresponding distances 


along two neighbouring geodesics in a congruence with uniu 
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tangent vector ]/ , Then is 'dragged' along by the 


congruence, that is 
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introducing an 


along V 
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wxtn 


orthonormal tetrad & 

X _ i/O. ^ 

^ - v we have 


parallelly transported 




vale re 


d b< 

A* 

ai S a 
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a K 

/V 
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IV t (j 

£, 6* 
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c . , & 
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(X a 


L cL 


1 / 1 / 


A solution of ( 4 ) will be called a jacobi field. There are 


clearly eight independent solutions. Since 1/ 


a. 


and 


Sv 


ft. 


are 


a geode 


solutions, the other six independent solutions of ( U ) may be 
taken orthogonal to l/°" . Then a is conjugate to p along 

esic / if, and only if, there is a Jacobi field along 
which vanishes at p and a . This may be shown as follows: 
the Jacobi fields which vanish at P may be regarded as 
generating neighbouring geodesics in the irrotational congruen 


-S 

I 


of all time-like geodesics through 


P • 


a ~r 


n ka at 0 0 0 "V 
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They may be written 


A. 


where 


sear 


P 


ff , 

Acs) 

d #f 


M fv 

el? ^ 

oL f) 


ik 

a 

M n 

2 . V 

f 7 

d s 


A 

A (s) 

will be 

positive 


P 


?n 


rxere v; 


00 


^ ^ \ /V\ - x \ 

a Jacobi i?iold vanishing at p and if* and only if fjl - 

t * i jf """St 


,o 


But f\ ($) 

A\ tv 

Therefore 

(9 


e^p(*C V els') 

Jp **•/» 




flL (^dt 

^ At A 


and 


d 


cl s' 


f) 

MA 


-a /? 

fw V h» 


Therefore 


<i Cfl) 

oLs ' 


is finite 


He nc e 


A 


is infinite where and only where 



r Q 


Thus the two definitions of conjugate points are equivalent. 
This also shows that singular points of congruences, are 
points where neighbouring geodesics intersect. 


p 
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2 


lor null geodesic congruences with parallelly transported 


tangent vector 


we may define the convergence 


o*, tZJ 


in Chapter p* i'his obeys 

fi c". ?* + 




id 




'^e define a singular point of a null geodesic congruence as 
one where P is infinite. 

ihe condition that the pressure is greater than minus 
one-third trie density may be stated more generally as 


condition (a) 


e->°> 


t > k T. 




ooserver w: 


4-velocity uJ , where 


- 


for any 


is the energy density in the rest-frame of the observer and 


a. 


is the rest-mass density. 




ondition (a) will be satisfied by a perfect fluid with 

/b 


it imolies 


8yVJo for 


cl exxL j 3 *t y ^ and 'ipiessxire p> - 5 p.. jl o 'v^ 

any time-like or null vector Y . Therefore by equations 

t 

(1) and (5) any time-like or null irrotational geodesic 


congruence mus 


t have a singular point on each geodesic within 


a finite affine distance. Obviously if the flow-lines form 
an irrotational geodesic congruence, there will be a physical 
singularity at the singular points of the congruence where 


the density and hence the curvature are infinite. This wni 
be the case if the universe is filled with non—rota.ting dust 


2,3 

7 


However, if the flow—lines are not geodesic (ie. non—vanishing 

pressure gradient) or are rotating, equation (ij cannot oe 
applied directly. 
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Spatially Homogeneous Anisotropic Universes 

'l'he Robertson-,valker models are spatially homogeneous 
and isotropic, that is, they have a six parameter group of 
motions transitive on a space-like surface, if we reduce the 


symmetry by considering models that are spatially homogeneous 
but anisotropic (thatr.is, they have a three parameter group 


of motions transitive on a space-like hypersurface) 
the matuer flow may have rotation, acceleration and 
Thus there would seem to be the possibility of non-s 


then 
shear. 
ingular 


4 

models. L. Shepley has investigated 
homogeneous model containing rotating 
there is always a singularity. Here 


one particular 
dust and has shown that 
a general result will be 


oroved 


There must be a singularity in every model which satisfies 


condition (a) and, 

(b) there exists a §y of motions 
universal covering space 

* bee section 5 


on the space or on 
which is transitive on at 


least one space-like surface but 


space-time is not stationary, 


(c) the energy—momentum tensor is 


that of a perfect fluid, 


• ' ^/ u r P) ^ ^ b * 

flow-lines and is uniquely defined as 
vector of the Ricci tensor. 



the 


-j— i—. 

\j CUi,g 


ent 



the time-like eigen- 
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ps?r 


T n 


R, the curvature scalar must be constant on a soace- 


like 


mu s o 


to // 


tne group 


Q 

herefore 


3 

surface of transitivity ft* of 
be in the direction of the unit time-like normal V 


Cu_ 


V 


CX. 




CX 


waere 




= f >0 






is an indicator - +1 if 


d. 


is east directed 


= -1 if fc . c_ is future directed 


mi 


'A 

V.V 


aen '^:b3 

J i' hu 3 14 congruence of geodesic xrrofuime~*lilCB 

H , V^V b >o 


vectors. By condition (a), 




herefore the congruence must have a singular point on each 


geodesic ( by equation 1) either in the future or in tne past, 
further, by the homogeneity, the distance along each geodesic 
rom bi to the singular point must be the same for each geodesic. 


-f 


thus if the surfaces of transitivity remain space-like 


, t hev 

7 */ 


must degenerate into, at the most, a 2-surface C which will 
be uniquely defined . Let M be the subset of the flow-lines 


of the matter which intersect 


c 


Let 


L be the non- 
iince there is a 


J 

empty subset of [i intersected by ri . 
group transitive on J4 ,4 must be (/ itself. Thus all the 


© Stephen Hawking, All Rights Reserved 

















Inches 
i i i 

t i i 

1 

1 , 


1 

i Lju 

A1 

cm 

i 

- 1 

2 


3 

U 

£ 


ilow-lines -through 
the density w 


must intersect the 2-surfh 


e O' 


Thu s 


infinite there an-.;, there will "he a 


physicai singularity. Alternatively if the surfaces of 
transitivity do not remain space-like, there mu 


be at 


U 


least one surface which is null 


call this S^. At 3^, = 0, 


Q 

K * 


f O ( if is zero, v/e can take any other scalar 

/ ^ 


polynomial in the curvature tensor and its covariant derivatives 
They cannot all he zero if space-time is not stationary;, he 

-^r 

introduce a geodesic irrotational null congruence on 3 with 
tangent vector L where £ R. * Then by equation 


<k /«x 

w -/ 
there will be a singular point of each null geodesic in '6 J 

within finite affine distance either in me future or in the 

past. The 2-surface of these singular points will be uniquely 

defined. The same argument used before shows that the donsiuy 

becomes infinite and there is a physical singularity. In tact 


mq C 

do hJ 


5 is a surface of homogeneity, the whole of S will be 


singulax* and it is not meaningful!to call il null or to 
distinguish this case from the case where the surfaces of 


m two ways, 


transitivity remain space-like. 

The conditions (a), (b), (c) may be w 
ondition (b) that there is a group o_ motions cnrougnout 

pace-time may be replaced by (b ) and (a). 


n 

A/ 
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/ 

(0 ) 


There is a space-like hypersurface 


H 


1 


in women t.ier 


"T"i Ci 


C\ 


are three independent vector fields A such that 

/VN 

I “3 °> • £ ft, , = o 

y*- bC ^ fccdLfc 

is one homogeneous space section. 


on 


tf 


That 


is 


there 


(d) There exist equations of state such that the Cauchy 

Li 3 

development of Ji is determinate. 

Then succeeding space-like surfaces of constant f\ 
are homogeneous and much the same proof can he given that there 

j 

are no non-singular models satisfying (a), (b), (c), (d). 

The only property of perfect fluids that has been used 
in the above proof is that they have well defined flow-lines 
intersection of which implies a physical singularity. Obvnousl 
however, this property will be possessed by a much more general 
class of fluids, for these, we define the flow vector as the 
time-like eigenvector (assumed unique) of the energy-momentum 
tensor. Then we can replace condition (c) on the nature of 
the matter by the much weaker condition (e). 


(e) If the model is singularity-free, the flow-lines 


form 


a Siii' 


ooth time-like congruence with no singular points with 


a line through each point of space-time. 


Condition (e) will be satisfied automatically if conditions 


(a) and (c) are. 





















'I'his proof rests strongly on the assumption of 

homogeneity which is clearly not satisfied by the physical ' 

universe locally though it may hold on a large enough scale. 

However it would seem to indicate that large scale effects 

like rotation cannot prevent the singularity. 

It is of interest to examine the nature ox the singularity 

in the Homogeneous anisotropic models since this is more 

likely to be representative of the general case than that of 

the isotropic models. It seems that in general the collapse 

S 

will be in one direction, that is, the universe will collapse 
down to a 2-surface. Near the singularity, the volume will be 
proportional to the time from the singularity irrespective of 
the precise nature of the matter. It also appears that the 
nature of the particle horizon is different. There will be 
a particle horizon in every direction except that in which 
the collapse is taking place. 


4. Singularities in Inhomop,’eneous Models 

6 

Lifshitz and Khalatnikov' claim to have proved that a 

da 

general solution of the field equations will not have a 
singularity, Their method is to contract a solution with a 
singularity which they claim is representative of the 
general solution with a singularity, and then show that it 
has one fewer arbitrary function than a fully general solution. 


f 
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Clearly their whole proof rests on whether their solution 
is fully representative and of that they give no proof. 
Indeed it would seem that it is not representative since it 
involves collapse in two directions to a 1-surface whereas 
in general one would expect collapse in one direction to a 
2 -surface. In fact their claim has "been proved false by 

7 


Penrose 


for the case of a collapsing star 


using tne norion 
of a 'closed trapped surface'. A similar method will be 
used to prove the occurrence of singularities in 'open ' 
universe models. 


5. 'Open' and 'Closed' Models 


The method used by Penrose to prove the occurrence 

of a physical singularity depends on the existence of a 

non-compact Cauchy surface. A Cauchy surface will be taken 

to mean a complete, connected space-like surface tnat 

intersects every time—like and null line once and once only 

Hot ail spaces possess a Cauchy surface; examples of those 

8 

that do not include the plane-wave metrics, the Godel model, 
and N.U.T. space. However none of these have any physical 
significance. Indeed it would seem reasonaole co demand. o_ 
any physically realistic model that it possess a Cauchy 
surface. If the Cauchy surface is compact, the model is 
commonly said to be 'closed'; if non—compact, it is said to 
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to be 'open' 


The surfaces, t = constant, in the Robertson- 


Walker solutions for normal matter are examples of Cauchy 
surfaces. If K = - 1 , they have negative curvature and it is 
frequently stated that they are non-compact. This is not 
necessarily so: there exist possible topologies for which 
they are compact. However, the following statements may be 
made about the topology of the surfaces t = constant. 

If the curvature is negative, It = -1, the universal 

X >1 A 

covering space is non-compact and is diffeomorphic to E . 


Any other topology can be obtained by identification of 

points. Thus any other topology will not be simply connected 

and, if compact, must have elements of infinite order in the 

fundamental group. Further if compact, they can have no 

12 

group of motions. 

If the curvature is zero, K = 0 , the universal covering 


o i x 

space is E . There are eighteen possible topologies. ^ 

compact they have a of motions and Betti numbers,- 3 ^ 


If 


—r 

= p 


b 2 = 5 . 


12 


If the curvature is positive, K = + 1 , the universal 

covering space is S . Thus all topologies are conroact. The 

12 

Betti numbers are all zero, 

Since a singularity in the universal covering space 
implies a singularity in the space covered, Penrose's method 
is applicable not only to spaces that have a non-compact bau 


■n n 
A-j- 












surface but also to spaces whose universal covering space 


has a non-compact Cauchy surface, 'thus it is applicable to 
nodeIs which, at the present time, are homogeneous and iso¬ 
tropic on a large scale with surfaces of approximate homo- 

JL. -SI— 

geneity which have negative or zero curvature. 


6 . The Closed Trapped surface 

,3 

Let t be a J-oall of coordinate radius r in a 3 -surface 

3 

(t = const.) in a Robertson-Walker metric with X = Q or - 1 . 
a 

Let q be uhe outward directed unit normal to T , the boundary 
3 3 0. 

of T , in H and let V be the past directed unit normal to 




Consider the outgoing family of null geodesics which 


T 2 orthogonally. At T 2 


oe: 


f 


f 


(K 


l 3 


, their convergence will 


A-; la 


UhXs * b + t tb ) 


where ^ are unit space-like vectors in orthogonal 


tj 3 


to q and to each other 


therefore 


f 


R 


K 


r /iXt 1 


If k7 O and K = 0 or -1, by taking r large enough, we may 

2 

make-p negative at T , Therefore, in the language of Penrose 


2 

ill*- A 


is a closed trapped surface. 


Another way of se.eing this is to consider the diagram 
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in which the flow—lines are drawn at their proper spatial 
distance from an observer. They all meet in the singularit 
at t = 0. If the past light cone of the observer is drawn 
on this diagram, it initially diverges from his world-line 


( p < a 


It reaches a maximum proper radius ( fi ~ 0 ) 


P < ). x\i reaches a maximum proper raaius ^ 

and then converges again to the singularity ( p > Q ). The 

5 

intersection of the converging light cone and the surface K 


gives a closed trapped surface T* 


2 


If the red-shift of the 


cuasi~steilar 3 >G 9 is cosmological then it will he beyond the 
point <9-0 if we are living in a Robertson-walker type 
universe with normal matter. However, the assumptions of 
homogeneity and isotropy in the large seem to hold out to the 
distance of $ 09 * Thus there is good reason to believe that 
our universe does in fact contain a closed trapped surface. 

It should he pointed out that the possession of a closed 
trapped surface is a large scale property that does not depend 
on the exact local metric. Thus a model that had local irregul¬ 
arities, rotation and shear but was similar on a large scale 
at the present time to a Robertson-Waider model.would have a 
closed trapped surface. 

Following Penrose it will he shown that space-time has 
a singularity if there is a closed trapped surface and : 

(f) E )/ 0 for any observer with velocity 

(g) there is a global time orientation 
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(h) the universal covering space has a non-compact Cauchy 

3 

surface H . 


PROOF 


Assume space-time is singularity free. 




JUJ0 L/ i? 


e o 


he 


set of points to the past of H'' that can be joined by a 

2 

smooth future directed time-like line to i' or its interior 

7 . 7 

T y . Let Ir be the boundary of F Local console rat ions show 

7 3 . 

that L - T is null where it is non-singular ana is 


generated by the outgoing family of past directed null geodesic 

2 

segments which nave future end-point on T 


oast ena-point 


where or before a singular point of the null geodesic 

2 - n y o 


congruence. Bince at 1*“ , the convergence, P 


ana 


Ck i fa 


since 


/? L*c~ >o 

0 ~k 


by (f), the convergence must 


become infinite within finite affine distance. 


"D 


will be compact being generated by a compact family of compac 

segments, hence will be compact. Penrose's method is 

y y . 

then as follows: approximate B arbitrarily closely by a 

3 3 

smooth space-like surface and project B onto H oy the 
normals to this surface. ,J -'his gives a many-one continuous 


its image B"' 


mapping of 3^ into Bince T? is compact, 

7 

must be compact. Let dlQ) be the number of points of B 
mapped to a point of H 5 . olj?) will change only at the 

3 

intersection of caustics of the normals with H. Moreover, 
oy continuity lit') can only change by an even number. 
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= 1 since this is the identity 


ana 


n 


This is a contradiction, thus the assumption that space-time 

is non-singular must be false. An alternative proceedure 

which avoids the slightly questionable step of approximating 

3 , 

B: oy a space—like surface is possible if we adopt condition 
(e) on the nature of the matter, then may be projected 


continuously one-to-one onto IV by the flow-lines. 


- * “Ti * 

-ms 


again beads to a contradiction since is compact and H 
is not. 

7 , 0 au c hy foriz on s 

uhe — above proof i t v;a . s •necessary — to d e mand tha t 


.3 


«4- 


3 3 

3 — o a —a-Aiauohy -surface- otnorv/is»-t.ne-whoAfG of B -might neb 


have been projcotoa ■■ onto-h 


3 


Wo will defin e— a -s emi — ^auohy - 


s urfaoo — (o ■ G. c )■■ ■ as a complot a connect e d space-liko surface 
w aam inters e cts evoiy- . time"like and null lin e -at most onoe . 


A hi 


/J j 

w-iJJl be a Cauchy ■ su rf a ce for — o i at -s - n o ar it 

t/ JL 


' W».ao -Arfery— i u • iill int e rsect every . t a. me--i.ikG ana . li^e 

thro~ag ' i *- - t- h»oe points -.— H owovop-, -- furth e r s w a y ■ thoro may b » 

. _ ok 

g egiono for w - haen it ic not a Cauohy ourfaoo. — Lot J -fe c o h e 

the boLHft&a - gy ■ of thoso po int c. — Cp' 1 , -if it oxiots; will be 


pall e d - tho Qauchy P-^r>i >7.nn tin w3^ 


muo u b a a 


s urf aoo .—^f urth e rmore if - condition - (f ) hold s the - null geod - o^a 
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■ mu s t -lyy i g i r . . ._4\f 


^ (thus—Q' J - is homeoraorph i e t o and i s 


£ompaot -» 
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If condition (i) — ho l ds ■ every null-geodesic gen e r • .to z- 


of (§ • jiao- at -loaot ono end-poin^ 


n 


hi - o mu□ u oo — in the 


directi - oa..away—from . H 


?- . since ■ tn - Ahe -di -root- i - on t ow ardc I - i 


3 


e - noh generator much he unbounded. — This houevor is inoosoidl e 


irr Q ^ " ° _ ghi , r; r Is a flfuichy sr 


8. Singularities in 'Closed* Universes 

y 'there is a singularity in every model which satisfies 


(a), (g) and (i). 


(i) 


0 3 


there exists a compact Cauchy surface H wnose uni 


normal I/ 0 " has positive expansion everywhere on * 

PS OOF 

For the proof it is necessary to establish a couple of 
lemmas. Assume that space-time is singularity—free. -ne 


; nr -?r*p L ; 


il;r b e 


following result is quoted without proof, it may 

derived from lemmas proved in reference 11. 

are conjugate points along a geodesic y 


If p ana q 
and X is a point on 


not in pc theny must have a conjugate 


point in 


Pp • 
























































An immediate corrollary is that if q is ohe iirst point 


along Y conjugate to p and y is in pq then y has no 
conjugate points in pq . Also since the result chat x 
has a conjugate point in pq can only depend on the values 
of in pq, any irrotational geodesic congruence including 

M a 

the geodesic ^ must have a singular point on ^ in pq. 

Thus if q is a point on Y? and f is the geodesic normal to 

Yi^ through q, then a point conjugate to q along Y cannot occur 

3 

until after a point conjugate to M . 

If Y? is a complete connected space-like surface which 

intersects every time—like and null line fr om a point p. 


v; e 


may define a function over as the square of the geodesic 
distance from p which is taken as positive if tue ^eoJ.e&j.'w 


he call 


time-like and negative if the geodesic is space-like, 
this the world function with respect to p. i*or uhe c-os^d 


et of values <5 ^ 


general multi-valued) function over 


df will be a continuous ( in 

5 


A time—like geouesic 


✓ from p will he said to be critical if it corresponds 

- * o Li * *) 


to a value of £> for which \ 

_ ju. ,3 

where are three independent vectors in n 

3 . , . . 

a critical geodesic must he orthogonal to II . a geodesic 


Clearly 
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1 12 


r 
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which is critical will be said to be maximal if it corresponds to a 


local maximum of 


Lemma 1. 


A geodesic 0 cannot be maximal for a smooth M if there is a 
point X conjugate to M but no point conjugate to q on Y in qp, 
where q is the intersection of X and M . 


and 


Let f and g be the Jacobi fields along / which vanish atX 
m m 

respectively. They may be written 


f 


n 

A(s) f/q 
mn 


n 


Then 

n 

any h 


(Y) ( 

L ( 


<//) 


g = B(s)g/q . 
m mn 

J.S. I 




i '■/ 



h 


l 

V - n 

since if it were negative for any h 


must be positive for 

by taking a = -b h h 

mb m n 


beyond q 


If it were zero 


beyond q, it would be possible to have a point y on 
conjugate to 'f- before a point conjugate to p 
X would be conjugate to P . This shows that the surface at q 
of constant geodesic distance from p lies nearer to f* in every 
direction than the surface o£ (j of constant geodesic distance from 
X does. Since X is conjugate to W the surface at q of 
constant geodesic distance from O lies closer to ^ in some 
direction than M does. Hence 




is not maximal. 
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4 > 


frVV 


g-r ^QQ 'ie’oi 't-h —p- — in the n i recti on K the — sv . r j .<<• > 

o f oonotant goodooio distance from p lie a oloscr to p thas r 


^ho — ourf qo O ' "14^- - do o o ■- 1 1 ho no fore — )f —m— non —max l m * 1 . 


If M 


5 


is compact or if the intersection of all time-like 

5 


and null lines with M" is compact, d* must have a maxima:.. 

3 

value, thus there must be a geodesic normal to K through p 


longer than 


, : e use this to prove another lenmu 


ox * 


Lemma 2 


If p lies to the future (oast) on a time-like geodesic 
v" through q, beyond a point z conjugate to Q,and there 
exists a compact Cauchy surface H' ; through q, tnen there 
must be another time-like geodesic from p to q longer than Y 
Let y be the last point conjugate to q on ^ bexore p• 
Let x be the nearest point to p conjugate to p in pq. Let 
r be a point in yx« Let be the set of point a wnikln n^ve 
a future (past) directed geodesic of length rq from c. 'fhen 
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cm 11 12 


F wi¬ 


pe a space—like hypersurface through r. 


Jjet F 




be the set of ooints which have at least one future (past 


\ 

J 


directed geodesic from q ox length greater than rq. then the 

4 % 3 4 

boundary ox F , d ci K . Since p is in F and since every 
past (future) directed time-like and null line from p intersects 


IT 5 


,A 


which is not in F , they must also intersect J . let 

5 -5 

be the inic-ersection of J and these lines. Since H is 


.3 


compact, L"’ must be compact. Consider the function CT k with 
respect to p over K . its maximum must lie in the compact 

“r 

region "L ?. But, by the previous lemma '( is not maximal, 

moreover, local considerations show that a singular ooint in 

-/ 

the surface Cw cannot be a maximum ox 0 ' . Thus the 

maximum value of (T must occur for a geodesic from p ortho- 
gonal to L . This must also be a goedesic from p to q of 
length greater than / . 

Using these two lemmas the theorem may be proved, oince 

g 

the future (past) directed normals to H are converging 

5 3 

everywhere on H , there must be a point conjugate to H a 

finite distance along eacn future (past) directed geodesic 


normal. 


Aj0 U 


be the maximum of these distances. Let p 
be a point on a future (past) directed geodesic normal at 


a distance greater than 


Consider the function 


VJ1 G hi 
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X 


respect to p over the compact surface Er 
geodesic from p normal to K at the point 
maximum, there must be a point conjugate 


Le 


oe the 


a, whereCThas its 
to ir alons A in co. 


But if there is no point conjugate to a along J\ in qp T then 
\ cannot be maximal by the first lemma. If however there 
is a point conjugate to q along / in qp, then there must be 
a longer geodesic from q to p by the second lemma, ihus 
is not the geodesic of maximum. length from H to p. ahis is 
a contradiction which shows that the original assumption that 
the space was non-singular must be false. 

i'his proof could also be used to show the occurrence 
of a singularity in a model with a non-compact Cauchy surface 
provided that the expansion of its normals was bounded away 


from zero and provided that the intersection of the 


o au c ny 


surface with all the time-like and null lines from a poin' 
was compact. 
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